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Abstract. The calculations of the p-process in the 0/Ne 
layers of Type II supernovae are quite successful in repro- 
ducung the solar system content of p-nuclides. They pre- 
dict, however, a significant underproduction of the rare 
odd-odd nuclide ^^^La. A model for the explosion of a 25 
Mq star with solar metallicity is used to suggest that Ve- 
captures on ^^^Ba may well be its most efficient production 
mechanism. The responsibility of an inadequate prediction 
of the ^'^^La and ^'^^La photodisintegration rates in the too 
low production of ^'^^La is also examined quantitatively. A 
detailed discussion of the theoretical uncertainties in these 
rates suggest that the required rate changes are probably 
f^ \ too high to be fully plausible. Their measurement would 
psj ■ be most welcome. They would help disentangling the rela- 
,-H \ tive contributions ofthermonuclear and neutrino processes 
, ^ to the -'^•^^La production. 

^^ ■ Key words: nuclear reactions, nucleosynthesis ~ solar 
system: general 
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1. Introduction 

The odd-odd neutron-deficient heavy nuclides ^^^La and 
isomeric ^^"Ta™ are among the rarest solar system 
species (no information exists for other locations), with 
i38La/i39La p» IQ-^ and isoTa^/i^^Ta « IQ-^. In spite 
of their scarcity, their origin has long been a puzzle. As 
initially claimed by Prantzos et al. (1990) and confirmed 
by Rayet et al. (1995; hereafter RAHPN), i*^OTa" appears 
to be a natural product of the p-process in the 0/Ne-rich 
layers of Type II supernovae (SNII). In contrast, ^"^^La is 
underproduced in all p-process calculations performed so 
far (e.g. Fig. 1 of Arnould et al. 2001). 

In view of the low ^^^La abundance, it has been at- 
tempted to explain its production by non-thermonuclear 
processes involving either stellar energetic particles (Au- 
douze 1970) or neutral current neutrino-induced transmu- 
tations (Woosley et al. 1990). The former mechanism is 
predicted not to be efficient enough, while the latter is 
found by Woosley et al. (1990) to be able to overproduce 
the solar ^'^^La/^'^^La ratio by a factor of about 50. This 
prediction has to be taken with some care, however, espe- 
cially in view of the qualitative nature of the evaluation. 

In a one-dimensional Z = Zq M = 25 Mq SNII 
model, ^'^^La is predicted to be produced only at peak 
temperatures around 2.4 ± 0.1 x 10^ K (Arnould et al. 
2001) from a subtle balance between its main produc- 
tion by ^'^^La (7 , n) ^^^La and its main destruction by 
^■^^La (7 , n) -'^^^La . The resulting abundances cannot ac- 
count for the solar system ^^^La amount. The same con- 
clusion holds for all the stars in the 13 < M/Mq < 25 ex- 
amined by RAHPN. This situation might of course just re- 
sult from inadequate astrophysics and/or nuclear physics 
inputs. On the astrophysics side, one might incriminate 
an uncertain prediction of the evolution of the thermody- 
namic conditions of the ^^^La producing layers during the 
explosion. Modifying these conditions is unlikely, however, 
to cure the ^'^^La underproduction. For any astrophysically 
plausible conditions, the SNII layers releasing the high- 
est ^^*La yields are also those overproducing even more 
significantly heavier p-nuclides as ^^^Dy, ^^^Er or ^®*Yb. 
Of course, it remains to be seen if the situation could be 



drastically modified if multi-dimensional effects were duly 
taken into account. On the nuclear physics side, one has 
to be aware of the fact that the ^"^^La yield predictions 
rely entirely on theoretical nuclear reaction rates. One is 
thus entitled to wonder about the sensitivity of the com- 
puted ^'^^La underproduction to the nuclear uncertainties 
that affect the production and destruction channels. One 
of our aims is to provide the first quantitative examination 
of this question (see also Arnould et al. 2001). 

Our second aim is to analyse on more quantitative 
grounds than Woosley et al. (1990) the possibility of pro- 
ducing ^"^^La at a level compatible with the solar system 
abundances through neutrino nucleosynthesis (i^-process) . 
To this end, the same M = 25 Mq SN model as consid- 
ered above is selected. The RAHPN p-process network is 
augmented with the neutrino and anti-neutrino charged- 
current inelastic scatterings off the network nuclei, and 
with the neutrino neutral current inelastic scatterings off 
the Ba, La and Ce isotopes. 

Section 2 analyses the purely thermonuclear p-process 
production of ^"^^La. A brief description of the adopted as- 
trophysical model and input nuclear physics is followed by 
the evaluation of the impact on the calculated ^"^^La yields 
of changes in the "''La (7 , n) "*La and "*La (7 , n) "'^La 
rates. We also examine the plausibility of the rate modifi- 
cations that are needed in order to avoid an unsatisfactory 
^^^La underproduction. The neutrino nucleosynthesis con- 
tribution to the p-nuclides, and in particular to ^'^^La, is 
discussed in Sect. 3. Conclusions are drawn in Sect. 4. 



2. -"^^^La and the p-process 

The p-process calculations reported here are based on a 
model for a Zq 8 Mq helium star (main sequence mass of 
about 25 Mq) already considered by RAHPN. It is evolved 
from the beginning of core He burning to SN explosion 
(for details, see Hashimoto 1995). We consider here as P- 
Process Layers (PPLs) 25 0/Ne-rich zones with explosion 
temperatures peaking in the (1.7-3.3) x 10^ K range. Their 
total mass is approximately 0.75 Mq. The deepest PPL is 
located at a mass of about 1.94 Mq, which is far enough 
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from the mass cut for all the nuclides produced in this 
region to be ejected during the explosion. 

The abundances of the s-process seeds for the pro- 
duction of the p- nuclides are calculated with the NACRE 
'adopted' ^s^e (a , n) ^^Mg rate (Angulo et al. 1999), as 
described by Costa et al. (2000; see the distribution la- 
belled i?i in their Fig. 2). 

The p-process reaction network includes some 2500 
stable and neutron-deficient nuclides with Z < 84. The 
n-, p- and a-capture reactions on all nuclei are consid- 
ered, as well as their (7,n), (7,p) and (7, a) photodis- 
integrations. The nuclear reaction rates are the NACRE 
'adopted' ones for charged particle captures by nuclei up 
to ^^Si (Angulo et al. 1999). For heavier targets, the rates 
predicted by the latest version of the Hauser-Feshbach 
code MOST (Goriely 2001) are adopted (the NACRE and 
MOST rates are available in the Brussels Nuclear Astro- 
physics Library ^ttp://www-astro. ulb.ac.be ) . This version 
benefits in particular from an improved description of the 
nuclear ground state properties derived from the micro- 
scopic Hartree-Fock method (Goriely et al. 2001), as well 
as from a more reliable nuclear level density prescription 
based on the microscopic statistical model (Demetriou & 
Goriely 2001). Finally, the experimentally-based neutron 
capture rates of Bao et al. (2000) are used. 

As in RAHPN, the computed abundance of the p- 
nuclide i produced in the PPLs of the Z© 25 M© star 
considered here is represented by its mean overproduc- 
tion factor (Fi) — (Xi) /Xi^Q, where Xi^Q is its solar mass 
fraction (Anders & Grevesse 1989), and 

'-- (1) 









X,,„_i)(M„-M„_i)/2, 



where Xi^n is the mass fraction of isotope i at the mass 
coordinate M„, Mp = ^ri>ii^ri — Afn-i) is the total 
mass of the PPLs, the sum running over all the PPLs 
(n = 1 corresponds to the bottom layer). An overproduc- 
tion factor averaged over all 35 p-nuclei is calculated as 
Fq = J2i i^i) /35, and is a measure of the global p-nuchde 
enrichment in the PPLs. So, if the computed abundances 
were exactly solar, (Fi) /Fq would be equal to unity for 
all i. Fig. 1 indicates that most of these factors lie in the 
0.3 to 3 range. The cases for which an underproduction is 
found are discussed by RAHPN or by Costa et al. (2000) 
for 92,94]vio and '^'^^^^Ru. Only the ^^^La problem is tackled 
here. 

We turn our special attention to the impact of al- 
terations in the MOST rates for ^^'''La (n , 7) ^■^^La and 
""^•^^La (n , 7) ^^^La . The Maxwellian-averaged cross sec- 
tions at T = 2.5 10^ K are predicted to be 123 and 
62 mbarn, respectively. These are used to calculate the 
reverse -'^^^La and ^^^La photodisintegrations of direct in- 
terest in the ^^^La production by the application of the 
detailed balance theorem. More specifically, we examine 
the extent to which the rate of -'^^^La (n , 7) ^'^^La has to 
be decreased and the one of ^'^^La (n , 7) -'^^^La increased in 
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Fig. 1. Normalised p- nuclide overproduction factors 
(Fi) /Fo {Fo = 42) obtained for the Z = Zq 25 M© model 
star with the nominal MOST rates and in absence of 
neutrino nucleosynthesis (open squares). For comparison, 
black squares designate the factors calculated by RAHPN 
with the use of a different set of nuclear reaction rates. The 
dotted horizontal lines delineate the 0.3 < (Fi) /Fq < 3 
range. 

order to bring the ^'^^La overproduction at levels compa- 
rable with those of the neighbouring p-nuclides. In the fol- 
lowing, the factors of decrease and increase will be noted 
^down and Flip, respectively. The corresponding reverse 
photodisintegration rates are decreased and increased by 
the same factors. In our numerical tests, Fdown and F^p 
are selected to vary independently in the 1 to 10 range. 

Fig. 2 shows the ratio i?i38 = {F(}^^Ce)) / (F("8La)) 
obtained in this test. The choice of ^'^^Ce as the normal- 
izing p-nuclide is dictated by the fact that it is produced 
in an amount close to the average value Fq, so that Riss 
gives a good representation of the ^■^'^La production by 
the p-process. It is seen that _Ri38 is pushed inside the 
0.3 < (F) /Fq < 3 range represented in Fig. 1 only for 



Fri, 



X F, 



20 to 25. It remains to be seen if such 



changes are physically plausible. 

Many of the input data necessary to calculate reac- 
tion rates with MOST have been measured for ^'^^La. This 
concerns in particular giant dipole energies and widths, 
or the neutron resonance spacings. In contrast, very lit- 
tle is known experimentally for ^^^La or ^^^La, so that 
their predicted neutron capture rates are likely to be less 
reliable. In order to evaluate the rate uncertainties, a se- 
ries of calculations have been performed with values of 
some basic input quantities differing from their nominal 
values to an extent which is considered as reasonable (in 
view of the trends and systematics for neighbouring nu- 
clei). Of course, the compatibility with the experimental 
constraints, if any, is always imposed. 

This analysis clearly demonstrates that the main 
sources of uncertainties lie in the nuclear level densities, 7- 
ray strength functions and neutron optical potentials. Re- 
lated errors in the capture rate predictions can amount to 
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Fig. 2. ^'^^Ce to ^^^La overproduction ratio, i?i38, f^or dif- 
ferent values of Fdown and F^p . 

about a factor of 2 for "'^^^La(n,7)"'^^^La at the temperature 
of relevance (2.4 10^ K) for the -'^^^La synthesis. In view 
of scarcer experimental information, larger uncertainties 
affect the ^^^La(n,7)^^^La rate. Even so, it is unlikely that 
the product i^down x ^up introduced above can realistically 
reach the requested value of about 20 to 25 in order for the 
underproduction of ^^^La with respect to ^'^^Ce not to ex- 
ceed a factor of 3 (see Fig. 2). Of course, only experimental 
determinations of the ^'^^La and ^^^La photodisintegration 
rates could really put this theoretical conclusion on a safe 
footing. Such experiments are feasible today (Vogt et al. 
2001) and are strongly encouraged. 

3. The production of -"^^^La by neutrino captures 

The calculations of Sect. 2 neglect the irradiation of the 
PPLs by the neutrinos streaming out of the proto-neutron 
star. In order to examine the impact of the neutrino inter- 
actions on the p- nuclide yields, and in particular on ^^^La, 
the reaction network used in Sect. 2 is augmented with all 
neutrino-induced reactions (charged and neutral currents) 
up to Kr (all rates are from Hoffman & Woosley 1992), 
as well as all charged-current neutrino and antineutrino 
scatterings up to Po. The neutral-current scatterings off 
nuclei heavier than Kr are considered only for the Ba, La 
and Ce isotopes. 

The neutrino and antineutrino scattering cross sections 
are averaged over supernova (anti)neutrino spectra that 
are approximated by zero-degeneracy Fermi-Dirac distri- 
butions corresponding to typical temperatures T^, = 8,5 
and 4 MeV for the Vx (where x stands for fj, and r 
(anti)neutrinos), De and Ve, respectively (Fuller & Meyer 
1995). The allowed transitions for the charged-current 
(anti)neutrino scatterings are treated fully microscopi- 
cally within the cQRPA approximation (Borzov et al. 
1995, Borzov & Goriely 2000) and on grounds of the 
ETFSI ground state description (Aboussir et al. 1995). 
The Gamow- Teller (GT) strength function accounts for 



the transitions to the daughter nucleus ground and low- 
lying states, to the GT resonance (GTR) and to the parti- 
cle continuum region above the GTR. The entire {N — Z) 
Fermi strength is contained in the isobaric analog reso- 
nance (lAR), its energy being taken from the experimen- 
tal systematics of Coulomb displacement energies. It has 
to be noted that the lAR and GTR usually appear as 
decay states for TV < Z nuclei. In such conditions, v- 
captures dominate i^-captures (the situation is the reverse 
in TV > Z nuclei, where the lAR and the main part of the 
GT strength are excited via z/-capture). 

As the average Ve{ve) energies are low [E <1Q McV), 
allowed transitions are expected to dominate the scatter- 
ing processes. For nuclei in the neighbourhood of ^"^^La, we 
just take into account the contribution from the forbidden 
transitions in a very rough way through the increase by 
25% of the microscopically calculated i^(i/)-capture cross 
sections for allowed transitions (Hektor et al. 2000). In 
particular, the resulting averaged I'e-capture cross section 
by i38Ba is estimated to be 7.5 lO"''^ cm^ at T^^ = 4 MeV. 
The neutral current I'-scattering contribution is also esti- 
mated in a rough way by assuming an approximate scal- 
ing of the (vxi^xf) cross section with mass number. At 
Ti,_^ — 8 MeV, the excitation cross section followed by 
neutron emission of the La isotopes is 1.2 10~^^ cm^. 

The energy- integrated number flux of neutrinos of type 
i^ at a radius rj (expressed in units of 100 km) is estimated 
in terms of the neutrino-sphere radius and temperature Ti, 
and of the i^-ncutrino luminosity L,j to be (Fuller & Meyer 
1995) 

*-)^^« ^ i°" ™";-<ion|^) O ^^ <^' 

In this expression, L^ is assumed to vary between about 
10^"'^ and 10^^ ergs s~^. This quite broad range of values 
is selected in order to accomodate to the many intrica- 
cies and uncertainties in the detailed neutrino fluxes ema- 
nating during the whole nucleosynthesis episode from the 
nascent neutron star underlying the 0/Ne-rich zones of 
interest here. Note that ^^ [Eq. (2)] depends on ry, which 
varies with time in all the considered layers. Let us also 
emphasize that no oscillation between neutrino species is 
taken into account. Would v^^s be converted to heavier 
species by matter-induced oscillations, larger i^(i?)-capture 
cross sections on bare nuclei could be expected (mainly 
due to the contribution of forbidden transitions) .R 

Fig. shows the impact of the neutrino interactions 
on the p-nuclide production for the two sets of typical lu- 
minosities (Fuller & Meyer 1995) i,.[10''^^ ergs s^^] = (3, 
4, 16) and (30, 40, 160) for [v^^Ve and Vx)- These two 
Li, combinations lead to increases of the ^"^^La produc- 
tion with respect to the case without neutrinos by factors 



^ This work was completed by the time the SNO collabora- 
tion has announced that their neutrino data were supporting 
the existence of neutrino oscillations (Giles 2001). A quanti- 
tative analysis of the impact of this exciting discovery on the 
//-synthesis of ^''^La has clearly to be postponed 
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Fig. 3. Comparison of the p-process abundances ob- 
tained with and without (squares) neutrino interactions 
(with standard MOST rates). Two sets of luminosities 
L^[10^^ ergs s'^] =(3, 4, 16) (circles) and (30, 40, 160) 
for (z/e,i/e, 1^2^) (triangles) are adopted for the calculations 
with neutrinos. 



of 4.8 and 36, respectively. This enhanced synthesis orig- 
inates entirely from the z/g-capture by ^'^^Ba, the neutral 
current scatterings on ^'^^La being found to have a negli- 
gible impact. Despite relatively similar cross sections, the 
^•^^Ba J/g-capture is found particularly efficient due to the 
large initial ^'^^Ba abundance (about 10 times the ^'^^La 
abundance). For the high luminosity set, the fg-captures 
also enhance the production of -'^^^In, ^-"^^Sn, ^^^Er and 
^^"Ta. Despite the numerous uncertainties still affecting 
supernova models and the neutrino physics in supernovae 
(spectra, luminosity, temperature, oscillation, interaction 
cross sections, ...), z^e-captures appear so far to be the 
most efficient production mechanism of the solar -'^'^^La. 

4. Conclusions 

Woosley et al. (1990) have considered as an 'intriguing 
possibifity' the production of ^'^^La by (charged current) 
z/e-captures on -'^^^Ba. We provide the first quantitative 
support to this possibility through the use of a 'reahstic' 
one-dimensional SNII model and of an extended network 
of nuclear and neutrino reactions the rates of which ben- 
efit from a series of improvements. Still, astrophysics and 
nuclear/neutrino physics uncertainties remain and forbid 
this conclusion to be as strong as one might want. 

The responsibility of an inadequate prediction of the 
^^^La (7 , n) ^^'''La and ^^^La (7 , n) ^^^La photodisintegra- 
tion rates in the too low production of ^^^La is also exam- 
ined quantitatively. Clearly, a suitable ^^^La production 
could be obtained by adequate changes in the nominal 
rates. However, it is concluded from a detailed study of 
the theoretical uncertainties in these rates that the level 
of the required changes is probably too high to be fully 
plausible. An experimental study of the above-mentioned 
key rates is no longer out of reach and is eagerly awaited. 



These measurements would certainly help greatly disen- 
tangling the relative thermonuclear and neutrino synthesis 
contributions to one of the rarest nuclides in nature. 
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